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Transitional pipe flow is modeled as a one-dimensional excitable and bistable medium. Models are pre- 
sented in two variables, turbulence intensity and mean shear, that evolve according to established properties of 
transitional turbulence. A continuous model captures the essence of the puff-slug transition as a change from ex- 
citability to bistability. A discrete model, that additionally incorporates turbulence locally as a chaotic repeller, 
reproduces almost all large-scale features of transitional pipe flow. In particular it captures metastable local- 
ized puffs, puff splitting, slugs, a continuous transition to sustained turbulence via spatiotemporal intermittency 
(directed percolation), and a subsequent increase in turbulence fraction towards uniform, featureless turbulence. 

PACS numbers: 47.27.Cn, 47.27.ed, 47.27.nf, 47.20.Ft 



The transition to turbulence in pipe flow has been the sub- 
ject of study for over 100 years 1 1 1, both because of its funda- 
mental role in fluid mechanics and because of the detrimental 
consequences of turbulent transition in many practical situ- 
ations. There are at least two features of the problem that 
make it fascinating, but also difficult to analyze. The first is 
that when turbulence appears, it appears abruptly (ll, and not 
through a sequence of transitions each increasing the dynam- 
ical complexity of the flow. Turbulence is triggered by finite- 
sized disturbances to linearly stable laminar flow |2-4|. This 
hysteretic, or subcritical, aspect of the problem limits the ap- 
plicability of linear and weakly nonlinear theories. The sec- 
ond complicating feature is the intermittent form turbulence 
takes in the transitional regime near the minimum Reynolds 
number (non-dimensional flow rate) for which turbulence is 
observed 1 1 , 5 - 7 |. In sufficiently long pipes, localized patches 
of turbulence (puffs) may persist for extremely long times be- 
fore abruptly reverting to laminar flow ISHTSl. In other cases, 
turbulent patches may spread by contaminating nearby lami- 
nar flow (puff splitting and slugs) 0111^16 -191. While min- 
imal models have been very useful in understanding generic 
features of intermittency in subcritical shear flows 120-24], 
such models do not capture the puffs, puff splitting, and slugs 
that are essential to the character of pipe flow. In this paper I 
argue that transitional pipe flow should be viewed in the con- 
text of excitable and bistable media. With this perspective I 
present models, based on the interaction between turbulence 
and the mean shear, that both capture and organize most large- 
scale features of transitional pipe flow. 

Figure [T] summarizes the three important dynamical 
regimes of transitional pipe flow. The left column shows re- 
sults from direct numerical simulations (DNS) (711251. Quan- 
tities are nondimensionalized by the pipe diameter D and the 
mean (bulk) velocity U. The Reynolds numbers is Re = 
DU /v, where v is kinematic viscosity. Flows are well repre- 
sented by two quantities, the turbulence intensity q and the ax- 
ial (stream wise) velocity sampled on the pipe axis. Specif- 
ically, q is the magnitude of transverse fluid velocity (scaled 
up by a factor of 6). The centerline velocity u is relative to 
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FIG. 1 . (Color online) Regimes of transitional pipe flow. Left column 
is from full DNS with 4x10^ degrees of freedom in a periodic pipe 
200D long. Flow is from left to right. Shown are instantaneous val- 
ues of turbulence intensity q and axial velocity u along the pipe axis, 
(a) Equilibrium puff at Re = 2000. (b) Puff splitting at Re = 2275. 
The downstream (right) puff split from the upstream one at an earlier 
time, (c) Expanding slug flow at Re — 3200. Right column shows 
corresponding states from the simple one-dimensional model ([3])-([6]) 
(d) R = 2000, (e) R = 2100, and (f) = 3200. 



the mean velocity and is a proxy for the state of the mean 
shear that conveniently lies between and 1. At low Re, as in 
Fig. [TJa), turbulence occurs in localized patches propagating 
downstream with nearly constant shape and speed. These are 
called equilibrium puffs |[2l[T6l[T3, a misnomer since at low 
Re puffs are only metastable and eventually revert to lami- 
nar flow, i.e. decay |8-15|. Asymptotically the flow will be 
laminar parabolic flow, {q = 0, = 1), throughout the pipe. 
For intermediate Re, as in Fig. [TJb), puff splitting frequently 



occurs [T6j [TtI [191. New puffs are spontaneously gener- 
ated downstream from existing ones and the resulting pairs 
move downstream with approximately fixed separation. Fur- 
ther splittings will occur and interactions will lead asymptoti- 
cally to a highly intermittent mixture of turbulent and laminar 
flow OITl. At yet higher Re, turbulence is no longer confined 
to localized patches, but spreads aggressively in so-called slug 
flow | 6, 17], as illustrated in Fig.[TJc). The asymptotic state is 
uniform, featureless turbulence throughout the pipe |7|. 

Models for these dynamics will be based on the following 
known physical features. At the upstream (left in Fig.[T]) edge 
of turbulent patches, laminar flow abruptly becomes turbulent. 
Energy from the laminar shear is rapidly converted into turbu- 
lent motion and this results in a rapid change to the mean shear 
profile |6, 26 |. In the case of puffs, the turbulent profile is not 
able to sustain turbulence and thus there is a reverse transi- 
tion |[6j|Z7 1 from turbulent to laminar flow on the downstream 
side of a puff. In the case of slugs, the turbulent shear profile 
can sustain turbulence indefinitely; there is no reverse transi- 
tion and slugs grow to arbitrary stream wise length 10 [13. On 
the downstream side of turbulent patches the mean shear pro- 
file recovers slowly 1271 , seen in the behavior of u in Fig. 1. 
Crucially, the degree of recovery dictates how susceptible the 
flow is to re-excitation into turbulence |26|. 

These are the characteristics of excitable and bistable me- 
dia 1281 |29l. In fact the puff in Fig.[TJa) bears a close resem- 
blance to an action potential in a nerve axon |30|. Linearly 
stable parabolic flow is the excitable rest state, turbulence is 
the excited state, and the mean shear is the recovery variable 
controlling the threshold for excitation. Thus, I propose to 
model pipe flow as a generic excitable and bistable medium 
incorporating the minimum requisite features of pipe turbu- 
lence. The models are expressed in variables q and u depend- 
ing on distance along the pipe. 

Consider first the continuous model 

qt^Uq:,=q{u^r-l-{r^ S){q - if) + fc, (1) 
Ut ^ Uux = ei{l - u) - e2uq - u^, (2) 

where r plays the role of Re. U accounts for downstream 
advection by the mean velocity, and is otherwise dynamically 
irrelevant since it can be removed by a change of reference 
frame. The model includes minimum derivatives, qxx and Ux, 
needed for turbulent regions to excite adjacent laminar ones 
and for left-right symmetry breaking. 

The core of the model is seen in the q-u phase plane in 
Fig. [2] The trajectories are organized by the nullclines: curve 
where u = and ^ = for the local dynamics (qxx = Qx = 
Ux = 0). For all r the nullclines intersect in a stable, but ex- 
citable, fixed point corresponding to laminar parabolic flow. 
The u dynamics with 62 > ei captures in the simplest way 
the behavior of the mean shear. In the absence of turbulence 
(q = 0), u relaxes to = 1 at rate ei, while in response to 
turbulence (q > 0), u decreases at a faster rate dominated by 
62. Values ei = 0.04 and 62 = 0.2 give reasonable agree- 
ment with pipe flow. (See the Appendix Sec. 1 1 a[ ) The q- 
nullcline consists ofq = (turbulence is not spontaneously 
generated from laminar flow) together with a parabolic curve 
whose nose varies with r, while maintaining a fixed intersec- 
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FIG. 2. (Color online) The distinction between puffs and slugs seen 
as the difference between excitability and bistablilty in Eqs. |lj-(|2|. 
Phase planes show nullclines at (a) r = 0.7 and (b) r = 1. The fixed 
point (1,0) corresponds to parabolic flow. In (b) the additional sta- 
ble fixed point corresponds to stable turbulence. Solution snapshots 
show (c) a puff at r = 0.7 and (d) a slug at r = 1. These solutions 
are plotted in the phase planes with arrows indicating increasing x. 



tion with g = at = 1 + (5, (^ = 0.1 is used here). The upper 
branch is attractive, while the lower branch is repelling and 
sets the nonlinear stability threshold for laminar flow. If lam- 
inar flow is perturbed beyond the threshold (which decreases 
with r like r~^),q\^ nonlinearly amplified and u decreases in 
response. 

The (excitable) puff regime occurs for r < Vc — £2/(^1 + 
62), Figs. [2|^a) and (c). The upstream side of a puff is a trig- 
ger front 1281 where abrupt laminar to turbulent transition 
takes place. However, turbulence cannot be maintained lo- 
cally following the drop in the mean shear. The system relam- 
inarizes (reverse transition) on the downstream side in a phase 
front 1 28 1 whose speed is set by the upstream front. Following 
relaminarization, u relaxes and laminar flow regains suscep- 
tibility to turbulent perturbations. The slug regime occurs for 
r > rc, Figs.[2|b) and (d). The nullclines intersect in addi- 
tional fixed points. The system is bistable and turbulence can 
be maintained indefinitely in the presence of modified shear. 
Both the upstream and downstream sides are trigger fronts, 
moving at different speeds, giving rise to an expansion of tur- 
bulence. A full analysis will be presented elsewhere. 

While Eqs. ([l])-([2]) capture the basic properties of puffs and 
slugs, the turbulence model is too simplistic to show puff 
decay and puff splitting. Evidence suggests that pipe turbu- 
lence is locally a chaotic repeller [31 1. Hence a more realistic 
model, within the two-variable excitability setting, is obtained 
by replacing the upper turbulent branches in Fig. [2] with a 
wedged-shaped region of transient chaos, illustrated in Fig. [3] 
Outside this region q decays monotonically. The model is 



F{q?+d{q^_,- 
< + ei(l-<) 



2q? + qT+i),<), 



(3) 

<-i), (4) 



where q'^ and u'^ denote values at spatial location i and time 
n. This model is essentially a discrete version of Eqs. ([l])-([2]). 
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FIG. 3. (Color online) Illustration of the discrete model, (a) Local 
dynamics in the u-q phase plane. Within a wedge-shaped region q 
undergoes transient chaos, while outside it decays monotonically to 
^ = 0. The region varies with R as indicated, (b) Map used to 
produce transient chaos. Parameter a (which depends on u and R), 
is the lower boundary separating monotonic and chaotic dynamics. 



except with chaotic q dynamics generated by the map F. 

The map F is based on models of chaotic repellers in shear 
flows 123^ 32 1 . Consider the tent map / given by 
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if q<Qi 
if Qi < g < 1 
if 1 < g < Q2 
if Q2 < ^ 



(5) 



with Qi = a/ (2 - 7) and Q2 = (4 + /3 - a - 7Qi)/(2 + P). 
Parameter a marks the lower boundary separating chaotic and 
monotonic dynamics, Fig.|3jb), while 7 sets the decay rate to 
the fixed point g = 0. For /3 > (/3 < 0) the map generates 
transient (persistent) chaos within the tent region. The map 
is incorporated into the pipe model by having the threshold a 
depend on u as well as on a control parameter via 



a = 2000(1 - O.Su)R- 



(6) 



The factor (1 — 0.8u) generates the desired wedged-shaped 
region, while 2000 sets the scale of R to that of Re. Finally, 
the map F is given by k iterates of /, i.e. F = f^; with k = 2 
used here. (See the Appendix Sec.| 1 b|) This has the effect of 



increasing the Lyapunov exponent within the chaotic region. 

The only important new parameter introduced in the dis- 
crete model is /3 since it quantifies a new effect - spontaneous 
decay of local turbulence for /3 > 0. Suitable values for others 
are: ei = 0.04 and 62 = 0.2 as before, 7 = 0.95, c = 0.45 and 
d = 0.15. (See the Appendix Sec. J_b ) As shown in Fig.[T] 
for /3 = 0.4 the model shows puffs, puff splitting, and slugs 
remarkably like those from full DNS. The model parameter R 
nearly corresponds to Reynolds number Re. 

While positive (3 is ultimately of interest, to better connect 
the two models consider first f3 negative, e.g. /3 = —0.4. A 
transition from puffs to slugs occurs as R increases and the 
wedge of chaos crosses the li-nullcline. One finds a noisy ver- 
sion of the continuous model in Fig. |2] (See the Appendix 
Fig. [8]) If splittings of turbulent patches occur, they are ex- 
ceedingly rare. At /3 ^ 0, (including even /3 = —0.1), chaotic 
fluctuations in q cause occasional splitting of expanding tur- 
bulence. Puffs at lower R are clearly metastable, persisting for 
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FIG. 4. (Color online) Three regimes of pipe flow from simulations 
of the discrete model ([3])-([6]). Space-time diagrams (time upward) 
illustrate (a) decaying puff at = 1900, (b) puff splitting at = 
2200, and (c) slug formation from an edge state at = 3000. For 
ease of comparison with published work on pipe flow, solutions are 
shown in a frame co-moving with structures. Turbulence intensity q 
is plotted with q = 1.8 in white. In (a) and (b) the scale is linear with 
^ = in black, while in (c) the scale is logarithmic with q < 10~^ in 
black. Dimension bars indicate space and time scales. The top space 
scale applies also to (b). 



long times before decaying. However, splitting and decay are 
unrealistically infrequent if f3 is too small. Setting /3 ^ 0.1 
gives realistic behavior, as seen in Fig. [T] where (3 = 0.4. 
(See also the Appendix Fig. [9]) Note that the splitting of ex- 
panding turbulent patches and the decay of localized puffs are 
caused by the same process - the collective escape from the 
chaotic region of a sufficiently large streamwise interval to 
bring about local relaminarization. This is precisely the sce- 
nario described by extreme fluctuations |33|. In the case of 
puffs, this results in puff decay, while in the case of splitting, 
laminar gaps open whose sizes are then set by the recovery of 
the slow u field. 

Figure |4] further illustrates how well the discrete model 
captures the three regimes of transitional pipe flow. Space- 
time plots show puff decay, puff splitting, and slug flow. In 
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FIG. 5. (Color online) Main figure is a bifurcation diagram for model 
turbulence in the thermodynamic limit. The turbulence fraction Ft is 
plotted throughout the transitional regime. The onset of sustained 
turbulence, via spatiotemporal intermittency, occurs continuously at 
Rc ^ 2046.2. Ft increases with R and saturates near R = 2800. 
Asymptotic regimes (laminar, intermittent, featureless) are labeled, 
along with corresponding transient dynamics (decaying puffs, puff 
splitting, slugs). The onset of featureless turbulence is not sharp as 
indicated by gray shading. Space-time plots illustrate the dynamics 
near the ends of the transitional regime (R = 2058, R = 2720, 
R = 2880) with q plotted in frames co-moving with structures (color 
map indicated, dark is laminar), (a) Mean lifetimes for decaying 
(circles) and splitting (squares) puffs crossing at i?x — 2040. (b) 
log-log plot of Ft versus R — Rc. Best fit to the filled (red) points 
determines Rc and the slope. 



Fig. Qa), a puff persist for only a finite time before abruptly 
decaying | 8 - 151. In Fig.Qb), puff splitting dominates the dy- 
namics |7, 16l|T71[T9l. New puffs are generated downstream 
from existing ones such that intermittent turbulent regions fill 
space. Compare especially with Refs. 7 and 19 Finally, in 
Fig. |4jc), a slug arises from a localized edge state (a low- 
amplitude state on the boundary separating initial conditions 
which evolve to turbulence from those which decay to laminar 
flow 1 18, 34-361). Compare especially with Ref.fTSl 

The remainder of the paper provides a global perspective 
of the transitional regime, obtained from extensive numerical 
simulation of Eqs. ([3}-(|4]) and summarized in Fig. [5] Tur- 
bulence fraction Ft serves as the order parameter, tracking 
the dynamics from the onset of intermittency through the ap- 
proach to uniform, featureless turbulence. A point is defined 
to be turbulent if q > O.ba and Ft is the mean fraction of 
turbulent points. 



There is a continuous transition to sustained turbulence, 
via spatiotemporal intermittency 1 1911241, at a critical value 
Rc ^ 2046.2. Below R^ the flow is asymptotically laminar 
and Ft = 0. Above Rc turbulence persists indefinitely and 
Ft > 0. This transition is associated with the crossing of mean 
lifetimes for puff decay and splitting shown in Fig.|5ja). Both 
decay and splitting are memoryless processes with exponen- 
tial survival distributions P ~ exp(— n/r(i?)), where r{R) is 
the i?-dependent mean time until decay or split. (See the Ap- 



pendix Sec. |2a| and Fig. [12]) The mean lifetimes vary approx- 
imately super-exponentially with R (131 [141 [Wi but neither 
is exactly of this form. Mean lifetimes cross at i?x — 2040. 
Above Rx 2in isolated puff is more likely to split than decay. 
As expected, even though individual turbulent domains may 
still decay, others may split, as seen in the spacetime plot at 
R = 2058. Due to correlations between splitting and decay 
events, the critical value Rc is not identical to Rx, but is very 
close to it (a difference of 0.3%). Fig. [sjb) shows that just 
above criticality. Ft ~ {R — Rc)^''^^, supporting that the tran- 
sition falls into the universality class of directed percolation 

EH. 

The ratio of turbulence to laminar flow increases through 
the intermittent region and at the upper end only small lam- 
inar flashes are seen within a turbulent background. Beyond 
R ^ 2800 laminar regions essentially disappear and Ft I. 
This occurs in pipe flow at Re 2600 |7|. The transition to 
featureless turbulence is not sharp, however, nor is the transi- 
tion from puff splitting to slugs. This upper transition will be 
addressed elsewhere, but the basic effect, common for bistable 
media, is already contained in Eqs. ([l])-([2l). For a range of r 
above the slug transition, (Tc < r ^ 0.91), turbulence does 
not expand to fill the domain in the presence of other slugs. 
Small laminar regions remain due to the recovery of the slow 
ii-field and this sets the scale for the laminar flashes at the 
upper end of the transition region in Fig. [5] 

I have sought to understand key elements of transitional 
pipe flow - puffs, puff splitting, and slugs - without appealing 
in detail to the underlying structures within shear turbulence. 
This approach is similar to that expounded by Pomeau 1201 , 
and considered elsewhere (23]. The important insight here is 
the close connection between subcritical shear flows and ex- 
citable systems. The view is that a great many features of 
intermittent pipe flow can be understood as a generic conse- 
quence of the transition from excitability to bistability where 
the turbulent branch is itself locally a chaotic repeller. I have 
introduced particular model equations to express these ideas 
in simple form. While phenomena have been demonstrated 
with specific parameters, the phenomena are robust. The chal- 
lenge for future work is to obtain more quantitatively accurate 
models, perhaps utilizing full simulations of pipe flow, since 
ultimately the fluid mechanics of shear turbulence (streaks and 
streamwise vortices) is important for the details of the pro- 
cess. More challenging is to extend this effort to other sub- 
critical shear flows, such a plane channel, plane Couette, and 
boundary-layer flows. These require non-trivial extensions of 
the current work because, unlike here, the mean shear profile 
cannot obviously be well captured by a simple scalar field. 
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Appendix: Supplemental Information 

1. Parameters Selection 

Here the parameter selection used in this study is discussed. 
No attempt has been made to determine precisely values such 
for the best fit to pipe flow. The models are not sufficiently 
quantitative that exact comparisons are called for at this time. 
Moreover, the phenomena presented in the paper are very ro- 
bust and for some parameters there simply is not a strong cri- 
terion to use to select precise values. The goal is to provide 
justification for the values used in the paper as well as insight 
into how the parameters control the dynamics of the models. 



a. Parameters for continuous model 

Only the two rates ei and 62 need to be determined. The 
value of 5 has little impact on the dynamics and has simply 
been fixed at 0.1. These parameters are determined by fitting 
to a typical puff from DNS as shown in Fig. [6] 
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FIG. 6. (Color online) Parameters a and 62 chosen to match a typical 
puff from DNS at Re = 2000. In the left plot £2 = 0.2 while a has 
values 0.02 (red), 0.04 (blue), and 0.06 (green). In the right plot 
ei = 0.04 and €2 has values 0.1 (red), 0.2 (blue), and 0.4 (green). 
r = 0.7. DNS is the irregular black curve. The unlabeled (blue) 
curves in the two plots correspond to the values of ei and €2 used in 
the paper. 

The left plot shows the spatial profile of model puffs for 
three values of ei, the parameter controlling the final relax- 
ation to parabolic flow. It is straightforward to select a reason- 
able value of ei from such a plot. Note, however, that a scaling 
of model length scale has been performed to plot model and 
DNS profiles on the same graph (model lengths have been 
multiplied by 0.225). This scaling of length is such that the 
sharp upstream edge of the model puff occurs over the same 
distance as in DNS. The upstream edge is largely set by 62. 
(If the scaling of space units between model and DNS where 



known for other reasons, then the spatial profile alone could 
be used to determine both ei and 62.) 

The right plot is used then to complete the determination. 
Here model puffs for different values of 62 are plotted in the 
u-q plane. The sharp upstream edge of a puff is the trajectory 
rising from parabolic flow atu = 1, q = and this is strongly 
affected by the value of 62. If 62 is too small then the trajectory 
is too step (u does not respond quickly enough). If 62 is too 
large, then q does not reach a sufficiently large value. 

The values of ei = 0.04 and 62 = 0.2 chosen for the sim- 
ulations presented in the paper were arrived at by varying the 
two values to get the best overall agreement in the spatial pro- 
file and phase portrait. 



b. Parameters for discrete model 

As stated in the paper, the two rates ei and 62 are taken 
to have the same values as in the continuous model. This is 
quite reasonable given the relationship between Eqs. (2) and 
(4). This leaves choosing the parameters k, (3, and 7 of the 
map F, and the parameters c and d. The role of each of these 
is discussed below. As each parameter is varied in the fol- 
lowing, the remaining parameters take the fixed values used 
in the paper: ei = 0.04, 62 = 0.2, k = 2, (3 = 0.4, 7 0.95, 
c 0.45, andd = 0.15. 

Parameter k: The parameter k effectively dictates how 
many iterates of the map / are used per time step of the model. 
The effect of the parameter k is shown in Fig. [7] where puff so- 
lutions are shown in the u-q phase plane for /c = 1 and k = 2. 
Model turbulence is more erratic for /c = 2 than k = 1. When 
k = 1, puff splitting and the transition to sustained turbulence 
occurs at a smaller value of R, but the fact that R is smaller on 
the top row of Fig.[7]only partially accounts for the difference 
between the top and bottom rows of Fig. [7] 

In addition to this visual comparison, there is the fact that 
the average slope for a unimodal map with stable chaos is lim- 
ited to A = 2 and this is an artificial constraint on the dy- 
namics that comes about from considering one-dimensional 
dynamics. (In the case of transient chaos the mean slope can 
exceed 2, but there is still a constraint relating the escape rate 
to the mean slope.) Taking k > lis equivalent to considering 
multimodal maps and removes the artificial constraint. 

Note that the model shows puffs, puff splitting and slugs 
even for /c = 1. These features are robust. However, the ad- 
ditional freedom in the chaotic dynamics by allowing k = 2 
permits the model to achieve a better representation of turbu- 
lent flow. I have not found that using k > 2 offers noticeable 
further improvement. 

Parameter (3: This parameter controls the leakage rate from 
the chaotic region of the map. Figure [S] shows examples of 
states for /3 = —0.4 and f3 = 0. With /3 sufficiently negative, 
as for /3 = —0.4 in the top row of Fig. [8] a transition from 
puffs to slugs occurs that is essentially just a noisy version of 
the continuous models shown in Fig. 2 of the paper. Note that 
the chaotic wedge first touches the ix-nullcline at = 1733 
and the transition from puffs to slugs occurs near to, but not 
exactly at, this value. If there are any splitting events they are 
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FIG. 7. Effect of parameter k. Top row /c = 1 and bottom row k — 2. 
In each case snapshots are plotted in the u-q phase plane. The left 
two plots show randomly chosen snapshots of solutions with three 
closely spaced puffs. The right-most plots are of solutions following 
a quench from high R and contain a large number of puffs. Top 
row is at = 1760 and the bottom row is at i? = 2000. Both 
values of R are close to the transition to sustained turbulence for the 
corresponding value of k. 



_i I I i_ 



_i I I i_ 



500 



500 



— 1 1 1 1— 



1 

:_ 


III. 


; 




f \ 



X 



500 



X 



500 



FIG. 8. (color online) Effect of parameter (3. Top row (3 — —0.4 and 
bottom row /3 = 0. In each case snapshots of solutions are shown 
on either side on the transition from localized puffs to expanding 
turbulence. Upper left: R = 1750. Upper right: R = 1780. Lower 
left: R = 1760. Lower right: R = 1800. 



very rare. 

With /3 ~ 0, as for /3 = in the bottom row of Fig. [8] the 
transition from puffs to slugs is mediated by splitting events. 
However, the splitting events are too rare for the model to re- 
alistically correspond to pipe flow. 

As illustrated in Fig.|9] with /3 ^ 0.1, model puff splitting 
is similar to pipe flow. There appears to be no strong basis to 
select any particular value of {3 based on a visual examination 
of the onset of splitting. The value /3 = 0.4 used in the paper 
is simply a representative value. 

Parameter 7: This parameter controls the monotone decay 
of turbulence q following exit from the chaotic region. Fig- 
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FIG. 9. (color online) Further effect of parameter /3. Left /3 = 
0.3 (R = 1980) and right /3 = 0.5 (R = 2180). In each case 
snapshots of solutions are shown just after a puff split. There is little 
to distinguish the cases. 




ure 10 shows puffs plotted in the u-q phase plane for two value 



FIG. 10. (Color online) Effect of parameter 7. Left shows 7 = 0.80, 
with R = 2400, middle shows 7 = 0.95, with R = 2000, and is the 
same plot as bottom middle of Fig. 2. Arrows indicate the relevant 
region of the u-q phase plane. The right-most plot is DNS of a puff 
at Re = 2000, exactly the same as in Fig. 1. 



of 7. For comparison the puff at Re = 2000 from Fig. [6] is 
repeated here. The smaller 7, the more quickly q decays and 
the less rounded the phase-space dynamics is in the lower left 
comer of the phase portraits. The best match of turbulent de- 
cay in the model will be at the largest possible value for 7. 
However, as 7 approaches 1 laminar flow becomes marginally 
stable to q perturbations and this is clearly unphysical. The 
value 7 = 0.95 was chosen as a compromise between the 
competing requirements of having 7 large but not too close to 
1 . An improvement could be likely be obtained by having 7 
be a function of u and R, but this introduces additional fitting 
parameters and is not done here. 

Parameters c and d\ These parameters are naturally thought 
of as arising from the discretization of the terms Ux and qxx 
in the continuous model. At present I am not aware of any 
compelling reasons to select c and d to particular values other 
than that they should be small {d must satisfy d^ 1 /2 for sta- 
bility reasons). The value of c was chosen to be less slightly 
less than 1/2. The value of d is then to be fixed. Based on the 
analogy with the continuous model, and the value chosen for 
c, it could be taken to be (i = . This is because discretizing 
the continuous model with a time step of At = 1 means that 
c = 1/ Ax, from an upwind discretization of the advection 
term. Then d will be l/(Ax)^ from discretization of the dif- 
fusion term. This would give d = 0.45^ = 0.2025. However, 
adjusting d downwards from this value to (i = 0.15 places the 
transition to sustained turbulence in the model at the critical 
Reynolds number for pipe flow. Having the model match this 
transition point seems preferable to setting it to the particular 
value 0.2025. Moreover, it is common to vary the diffusive 
coupling constant in studies of coupled-map lattices. 
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FIG. 1 1 . (Color online) Effect of parameters c and d. In each case 
snapshots of solutions are shown just after a puff splitting. Splitting 
is not very sensitive to c and d around values used in the paper. Upper 
left: c = 0.45, d = 0.1, R = 2480. Upper right: c = 0.45, d = 
0.2025, R = 1960. Lower left: c = 0.35, d = 0.15, R = 2120. 
Lower right: c = 0.55, d = 0.15, = 2080. 
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FIG. 12. Lifetime statistics for decaying and splitting model 
puffs, (a) Time series of total energy E for three puffs illustrat- 
ing abrupt decay at unpredictable times {R = 1800). Bottom row 
shows exponential (memoryless) probabilities P for (b) decaying 
puffs {R = 1800, 1850, 1900, 1950, 2000) and (c) splitting puffs 
{R = 2060, 2090, 2120, 2150, 2180). 



To emphasize that qualitative features of the model do not 
strongly dependent on the parameters c and d (as long as they 



are reasonably small), Fig. [TT] shows some splitting puffs for 
different values of c and d. The reason for focusing on split- 
ting puffs is these best show the fidelity of the model. Puffs 
and slugs are easily obtained. In each case a puff was gen- 
erated and R was increased slowly until a splitting occurred. 



The value of R are given in the caption. 

2. Details of Numerical Study in Figure 4 



a. Decay and Splitting Statistics 



Figure 
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a) shows time series of total energy E = ^ • qi 
for puffs 'mm three initial conditions in which q is varied 
at one space point by less than 10~^. While the dynam- 
ics are deterministic, abrupt decay occurs at unpredictable 
times. From a large number of such simulations lifetime 
statistics can be generated. Specifically, simulations were 
initiated from small perturbations designed to trigger puffs. 
A random number generator was used to introduce a small 
amount of randomness into each perturbation. Simulations 
were then run for 1000 time steps to allow puffs to equili- 
brate. If puffs had not decayed, the puffs were used as initial 
conditions for simulations for decay statistics. Figure [T2jb) 
shows representative survival probabilities P of a puff lasting 
at least time n. Each distribution corresponds to 4000 realiza- 
tions. The survival functions are exponential (memoryless), 
P rsj exp(— n/r(i?)), where r{R) is the i?-dependent mean 
lifetime until decay. 

Puff splitting is similar. Initial conditions were generated in 
the same way only here equilibration simulations were run for 
1400 time steps because at the largest R in the study 1000 time 
steps was not quite enough to remove all equilibration effects. 
A puff was defined to have split once two turbulent peaks are 
separated by least 80 grid points. Figure [T2jc) shows repre- 
sentative survival probabilities P for a puff to last at least time 
n without splitting. The distributions are again exponential, 
P ~ exp(— n/T(i?)), showing that model splitting is indeed 
memoryless with a mean splitting time r{R). These lifetimes 
are plotted in Fig. 4(a) of the paper. 



b. Turbulence Fraction 

The turbulence fraction Ef serves as the order parameter for 
the onset of sustained turbulence. A point is defined to be tur- 
bulent if q > hia where k, sets a threshold relative to the lower 
boundary separating chaotic and monotonic decay of the map 
F. Et is the mean fraction of grid points in the turbulent state. 
Means have been computed from 4 independent simulations. 
For R near the critical point, simulations of 8 x 10^ time steps 
were run on grids of 12 x 10^ points. The standard deviation 
from the 4 independent simulations is comparable to the point 
size in Fig. 4(b) of the paper. For R > 2200 simulations for 
10^ time steps on grids of 10^ points were more than suffi- 
cient. In Fig. 4 of the paper. Ft is plotted for = 0.5. While 
the exact value of Ft depends on threshold hz, the extent of the 
intermittent region and critical scaling at onset do not. 
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